Monte-Carlo simulations of the violation of the fluctuation-dissipation theorem in 

domain growth processes 
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Numerical simulations of various domain growth systems are reported, in order to compute the 
parameter describing the violation of fluctuation dissipation theorem (FDT) in aging phenomena. 
We compute two-times correlation and response functions and find that, as expected from the exact 
solution of a certain mean- field model (equivalent to the O(N) model in three dimensions, in the limit 
of N going to infinity) , this parameter is equal to one (no violation of FDT) in the quasi-equilibrium 
regime (short separation of times), and zero in the aging regime. 



The study of aging phenomena is currently the subject 
of many efforts, since this kind of behaviour, for which a 
given system remains out of equilibrium at all available 
times, is present in many systems of interest, like spin 
glasses or structural glasses 0. When concerned with 
the dynamics of a given system, it is usual to study the 
correlation function of an observable A, 



C(t,t') = (A(t)A(t')), 



(0.1) 



((.) denotes an average over thermal noise) and the con- 
jugated response function 



R{t,t') 



, 9A{t) 
'dh{t') 



(0.2) 



where h is an external field applied at time t 1 . Then, 
at equilibrium, these two-times quantities satisfy time 
translational invariance (TTI: the functions depend only 
on the difference of the two times t — t 1 )) and the fluctua- 
tion dissipation theorem (FDT) relating correlation and 



response by R(t — t' 



1 dC(t 
T 



-t') 



ill' 



On the other hand, 



for aging phenomena, since the dynamics is out of equi- 
librium, such equilibrium properties are not expected to 
hold. In the context of mean-field spin glasses, Cuglian- 
dolo and Kurchan have proposed the general following 
scenario, in the limit where the times t and if go to in- 
finity [Q: for small time differences ((t — if) ft! << 1), 
the system is in quasi-equilibrium, and the equilibrium 
properties hold; however, if t — t' is not small with respect 
to t', the study of two-time quantities reveals that it is 
not at equilibrium (C(t, t') depend explicitly on t and t'). 
Moreover, they have proposed to measure the violation 
of FDT by the function X(t, if) where 



R{t,t') 



X(t,t') dC(t,t') 



T 



dt' 



(0.3) 



with the important assumption (afterwards supported by 
the study of many different cases, see for example (| [7J) 
that, as t and t' go to infinity, it becomes a function of 
times only through C(t,t'): 



R(t,t') 



X{C) dC{t,t') 
~f dt' ' 



(0.4) 



This X(C) has moreover received an interpretation in 
terms of effective temperature [gj . In the high tempera- 
ture phase of any system, X is equal to 1 since the system 
equilibrates and the equilibrium properties hold. In the 
low temperature phase where aging phenomena appears, 
violations of FDT can be quantified by its departure from 
1. In simulations or experiments, it is more convenient 
to look at an integrated response function: the system 
can be quenched under a magnetic field, which is cut 
off after a waiting time t w (the relaxation of the magne- 
tization is then measured, and found to depend on the 
waiting time), or it is quenched under zero- field, and a 
field is applied after t w . In this second case, the growth 
of the zero-field-cooled magnetization 



M(t + t w ,t w ) 



R(t + t w ,s)h{s)ds (0.5) 



tw 



is observed. The quasi-FDT relation (0.3) allows then to 
write (for a constant field) 



-M(t + t v 
h 



t+t„ 



X(t + t w ,s) 



dC(t + t w ,s) 



ds 



ds 



which, in the limit of large t w , gives 

' tw ? t W ) 



T -M { t 



Jc 



C(t+t w ,t w ) 



X(C)dC. 



(0.6) 



(0.7) 



Then 

T 



if FDT is satisfied, we obtain a linear relation 
j^M(t + t w ,t w ) = 1 — C(t + t w , t w ), independently of the 
system, while a deviation from this straight line in a M 
versus C plot indicates violation of FDT and gives infor- 
mations on X: different systems can have different types 
of violation of FDT. This kind of M-versus-C plot has 
been used to compute the value of X in the aging regime, 
analytically for various mean- field models [||J|, and us- 
ing numerical simulations for the mean-field Sherrington- 
Kirkpatrick model ||, for the 3-dimensional Edwards- 
Anderson model 0] (a more realistic spin glass), for the 
p-spin in finite dimensions While, for the j? = 2 spher- 
ical p-spin model, equivalent to the O(N) ferromagnetic 
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model in three dimensions, X is zero it is found to 
be constant for p > 3, and a non-trivial function of C for 
the Sherrington-Kirkpatrick and the three-dimensional 
Edwards- Anderson model. An numerical investigation of 
a glass forming binary mixture (in three dimensions) has 
also been made recently , with the result of a constant 
value of X. 

In this letter, we report numerical simulations of var- 
ious domain-growth systems (for a review on such sys- 
tems, see [|l0|), for which it is expected || that X is 
zero in the aging regime. We examine Ising ferromag- 
netic systems in two and three dimensions at various 
temperatures, and with conserved or non-conserved order 
parameter. We also make a simulation of the Edwards- 
Anderson model in three dimensions, to show the striking 
difference of behaviour. 

We consider Ising spins s.; on a square or cubic lattice 
of linear size L, with ferromagnetic interactions. Start- 
ing from a random configuration, we quench the system 
at time to temperature T and let it evolve according 
to Glauber dynamics, with a single-spin-flip algorithm 
(we will also consider later soft-spins evolving through 
a Langevin equation). We then measure the spin-spin 
correlation function 



1 N 



(0.8) 



for a unperturbed system. It is known that this correla- 
tion function exhibits two time regimes: for t — t' << t' 
(for simplicity we take t' < i), it decays rapidly from 
1 = C(t',t') to qEA = ™ 2 , m being the magnetization at 
temperature T; then, for more separated times, it scales 
like L(t)/L(t'), where L(t) is the characteristic size of the 
domains at time t. We also check that the domain sizes 
remain much less than L, thus ensuring that finite size 
effects are not significant. At a certain waiting time t w , 
we take a copy of the system, to which a small, constant 
magnetic field is applied. We then measure the staggered 
magnetization 



M{t + t w ,t % 



1 



N 
i=l 



t)hi). 



(0.9) 



For spin glasses, the applied field can equivalently be 
taken uniform or random, since the interactions between 
spins are random. Taking an uniform field allows to avoid 
averaging over the realizations of the field. On the other 
hand, for a ferromagnetic system, the action of a uniform 
field is to favor one of the phases, which will grow faster. 
The correct quantity to measure is therefore the response 
to a random field: the staggered magnetization (0.S) [jj. 



For simplicity, the random hi are taken from a bimodal 
distribution [hi — ±h). The staggered magnetization is 
averaged over the realizations of hi , and we checked linear 
response using various values of h (typically from 0.01 to 
0.2). The sizes used are L — 600 in two dimensions, and 
L = 80 in three dimensions. 

To compare the various curves, obtained for various 
systems, temperatures and waiting times t w , we look at 
the plots of TM(t + t w , t w )/h versus C(t w + t, t w ). We 
first made some runs at high T: in this case, the system 
reaches quickly equilibrium, with TTI (C{t w + t,t w ) = 
C eq (t), M(t + t w ,t w ) = M eq (t)) and we checked that 
FDT holds (TM eq (T)/h = l~C eq {t)). For temperatures 
below the transition temperature, a dependence on t w 
appears in C and M (violation of TTI), corresponding 
to the growth of domains of the two competing phases. 
We observe as expected two time regimes PI: 

• for times t smaller than t w , the two-times quantities 
do not depend on t w , and FDT also holds: TM(t + 
t w ,t w )/h = 1 — C(t w +t, t w ). This happens at large 
values of C (close to 1) and small values of M. 

• for larger times separation, we observe aging in the 
correlation function, and also clearly a deviation 
from FDT. 

We show the data in figure (1), (2) and (3) for the var- 
ious systems, and for various waiting times. In the aging 
part, we see that the M versus C curves are in fact get- 
ting flat, except at small t w . A closer look at the data for 
the aging part shows that: (i) for larger t w , the plateau 
reached by the magnetization is lower, and (ii) for a fixed 
t w , the magnetization first grows (like 1 — C(t w +t,t w ), 
this is the non-aging part), then saturates, and eventu- 
ally goes slowly down again, this last effect becoming less 
important as t w grows, with a flatening of the curves (the 
slope of this part of the curves decreases as t w increases). 
We can explain these effects in the following way: after 
t w , the domains have reached a certain typical size, and 
the domain walls have a certain total length. The effect 
of the random field is then to try to flip some spins; this 
flipping will be easier at the domain walls, since the spins 
there are less constrained by their neighbours. Therefore 
we have two contributions to the staggered magnetiza- 
tion: one from the bulk, and one from the domain walls. 



*in two dimensions, a random field destroys the long range 



order (see jll] for a review on the Random Field Ising Model) ; 
however, the instability destroying it appears only for domain 
sizes growing exponentially with 1/h [pi, so that this effect 
is not important as long as we work with small enough fields 
and at times not to long 

'we stress that we are interested in long time limits, since 
the X{C) function is defined as such; nevertheless, we already 
can observe two distinct regimes with finite times, and deduce 
the limit of interest. 
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As time evolves, the domains grow and the total length 
(or surface, in three dimensions) of the domain walls de- 
creases. Therefore, the contribution from the interfaces 
decreases. On the other hand, the contribution of the 
bulk will be rather independent of t w , since the effect on 
a random field on a domain of + spins or on a domain of 
— is the same on average. The total staggered magneti- 
zation is thus decreasing when t w increases, and also, at 
t w fixed, as t grows (after the initial growth, when the 
field is switched on). In the limit of large t w , the effect 
of the bulk becomes relatively more important, and we 
observe the flatening [J 
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FIG. 1. TM(t + t w ,t w )/h versus C(t + t w ,t w ) for two di- 
mensional domain growth (T c — 2.27), at temperatures (from 
top to bottom) T = 1.7 and t w = 200, 400, 800, 2000, 
T = 1.3 and t w = 800, T = 1 and t w = 800. The straight line 
is M — 1 — C: we see that FDT holds at short times t, and 
the violation of FDT with X = at longer time separation. 

Note : the reciprocity relations, which state that, for 
two observables A and B, the correlations CUs(i,i') = 
(A(t)B{t')) and Cba^-,^) are equal, are also an equilib- 
rium theorem, and therefore are not expected to hold 
for aging dynamics. For a field <f> evolving according 
to a Langevin equation, where the force at time t is 
F(t), it can be shown jlij that, even if the asymme- 
try A(t, if) = (F(t)<f>(f)-F(f)<f>(t)) goes to zero for long 
times, the integral L A(t, t')dt' has a finite limit as t goes 
to infinity, if the system is out of equilibrium. Following 
a suggestion by S. Franz, and slightly modifying the sim- 
ulation program, we checked that this fact, derived us- 
ing the Langevin equation, also holds for a Monte-Carlo 
dynamics, where the field is replaced by the spins, and 
the role of the force is played by the local field acting 



J we have checked by a direct visualization of the spins that 
this is indeed what happens: at short times, the majority of 
the spins flipped by the random field are on the domain walls, 
this fraction going then down as the domains grow; we will 
also see that this effect due to the motion of domain walls is 
not present for the Edwards- Anderson spin-glass. 



on the spins. We therefore mention this integrated quan- 
tity, which could also be of interest in the studies of aging 
phenomena. 
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Same as figure 
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(1) for non conserved 



parameter in three dimensions, T = 2.5 (T c ~ 3.5), 
t w = 100, 300, 600, 1000, 1500. 
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FIG. 3. Same as figure (1) for conserved order parame- 
ter in two dimensions, T = 0.8 and from top to bottom 
tw = 100, 200, 400, 600, 800, and in three dimensions (lower 
symbols), T = 2, t w = 100, 200, 300, 400. 

Langevin equation : since similar results were ob- 
tained independently by C. Castellano and M. Sell- 
itto |f| for a system of soft-spins evolving through a 
Langevin equation, we also mention briefly this case, and 
show in figure (4) an example of the results that can be 
obtained with a system of this type: we simulate soft- 
spins on a square lattice, with a quartic potential confin- 
ing them to the vicinity of its minima +1 and — 1, and 
evolving through the discretized Langevin equation 

s(i,j,t + l) = s(i,j,t) + {s(i + l,j,t) + s{i - l,j,t) 

+ s(i,j + l,t) + s(i,j - l,t) - 4* s(i,j,t) 
+ s(i,j,t)-s(i,j,t) s )*h + r)(i,j,t) , (0.10) 

where s(i,j,t) is the value of the spin at the lattice site 
(i,j) at time t, r\ is a gaussian noise with zero mean and 
variance 2Th, h being the used time-step. We proceed 
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by parallel updating of the field, and, at t = 0, the s(i,j) 
are taken as independent random variables uniformly dis- 
tributed between —1 and 1. Again, at t w a random field 
is switched on and the staggered magnetization and the 
correlation are measured. 
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FIG. 4. Same as figure (1) for soft s pins on a 2-dimensional 
square lattice, evolving through (0.1C), with, from top to bot- 
tom, T = 1 and t w = 200, 600, T = 0.33 and t w = 200, 600. 

All these simulations clearly show that the parameter 
X is zero for these domain-growth systems. This flaten- 
ing of the integrated response shows that the long-term 
memory of such systems is in fact weak ||: the aging 
phenomena is essentially in the correlations, while it is 
also important for the response in spin-glasses. 

In figure (5), we indeed show the obtained data for 
an Edwards Anderson system in three dimensions, with 
Hamiltonian 



H 



Jij s i s j 

m 



(0.11) 



where the sum is over nearest neighbours, the spins Si are 
Ising spins, and the couplings Jy are quenched random 
variables, taking values +1 or —1 with equal probability. 
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FIG. 5. Same as in figure (1) for the Edwards Anderson 
model, T = 0.7, from top to bottom t w = 1000, 600, 300, 100, 
compared to the data for the domain growth in two dimen- 
sions at T = 1.7, t w — 1000 (highest plateau), and in three 
dimensions at T = 2.5, , t w — 1000 (lowest plateau). 



We simulated a system of linear size L = 80 at T — 0.7. 
Although no precise conclusion can be drawn as to the 
form of the function X(C), since the obtained curves still 
show a dependence on t w , it is quite clear (as was shown 
in Q) that they tend to a certain non-trivial curve, very 
different from the case of domain growth systems, like 
the comparison of figure (5) shows. Let us remark that 
curves similar to the ones obtained for the EA spin glass 
have also been obtained for the p-spin model in three 
dimensions in || an d for the mean-field version of (3.11), 
the Sherrington-Kirkpatrick model j(| . 

To conclude, we have reported measurements of the 
violation of the fluctuation dissipation theorem in some 
systems exhibiting domain-growth, and found that, as 
expected but shown only in one particular case, the pa- 
rameter X describing it is equal to zero in the aging phase 
(and of course to 1 in the quasi-equilibrium regime, where 
FDT holds). In the interpretation of J8|, this means that 
the effective temperature is the temperature of the heat- 
bath in the quasi-equilibrium regime (corresponding to 
the fast relaxation of the spins in the bulk of the do- 
mains), while it is infinite in the coarsening regime, which 
corresponds to the dynamics of the domains themselves 
(see S, paragraph IV-C for a detailed discussion). It 
should also be noted that this behaviour shows a ten- 
dency of the long-term memory to disappear, in contrast 
with spin glasses or glasses. 
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